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Determinant Calculations Using Random Walk Worldline Loops 
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We use statistical ensembles of worldline loops generated by random walk on hypercubic lattices to calculate 
matter determinants in background Yang-Mills fields. 



Introduction The World Line Formalism 
(|^ - see also Q and references therein) permits 
a transparent discussion of matter determinants 
and the related 1-loop effective action in some 
background. We here consider a lattice imple- 
mentation of this method. Consider the bosonic 
determinant in a gauge field background 
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The (Euclidean) worldline sum above is over 
all closed paths (trp indicates path ordering). 
We use Random Walk (RW) paths on a D- 
dimensional hypercubic lattice Ajj Q and dis- 
cretize T using At = T/L a'^/2D, then: 



C{U,x) = Y,) 
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L ^ ^ Wl{x)} 1Gwl(x) 

DimensionfuU quantities are understood as given 
in lattice units a. {lljl{x)} is the set of all closed 
lattice paths with length = L obtained by 
RW starting and ending at x. The free RW "mea- 
sure" (21?) is implemented by the actual proce- 
dure, the mass dependence e~™ is kept explicit. 
Ui are link variables. For detail see [Q. 
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The expansion according to the loop lengths 
(|^) is equivalent to the usual hopping parameter 
expansion for the logarithm of the determinant 
(see, e.g., ||]), e.g., in the bosonic case: 

= (to2 + 2D)-^ ~ (2D)-^e-™'* 

The number of different loops of length L in- 
creases exponentially with L, therefore one usu- 
ally stops the expansion at a rather low order. 
Since loops L are typically relevant at scale 
we then may completely miss physical effects at 
large scales (or need to use rough lattices). Our 
procedure amounts instead to sampling loops sta- 
tistically in a wide interval of lengths (see also 
0). By construction, the loops appear with the 
correct probabilities. We thus trade the system- 
atic errors for statistical ones, hoping to repro- 
duce also effects at "large" scale, albeit within 
statistical uncertainties. For a quantum gauge 
,^ field vacuum the (quenched) vacuum energy ob- 

^ ' tains by averaging (|l|) with the Yang-Mills action. 

The lattice loop ensemble The fac- 
^ tr Y\ Ui to'' C^Dy^ in (0) is reproduced by RW and 

we only need to weight the loop contributions 
tr (Jljg^^ Ui) with L~^e~"^ 20. We specify a 
maximal number of trials A/l for each given num- 
ber of steps L. We first generate an {w(a;^ = 1)} 
ensemble of loops by starting the random walk 
at the point a;^ = 1 of the periodic lattice and 
collecting the loops which after L steps have re- 



turned to = 1, independently on self-crossings 
or retracing. Using this {uj{Xfi — 1)} ensemble we 
define £{A,1). If we want to calculate £{A,x) 
we shift the {uj{Xf^ = 1)} ensemble by 
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whole. The fuU loop ensemble - the x-summation 
J2x£A (0) ^ approximated by considering 
sufficiently many random translations of each of 
the loops of the {w(a;^ = 1)} ensemble (this re- 
produces the correct multiplicities in the k - ex- 
pansion [^). We have for the loop frequencies: 

_ 2nL ^ T-D/2 (r,^ 

i2n/D)D/^JVL L»l ■ ^ ' 

with TiL the number of loops generated at L. The 
factor 2 appears because L is even. The second 
equality in (2) is the well known result for un- 
biased RW |2]. Since the number of loops de- 
creases with L, if large loops are important we 
may have large statistical errors. One improves 
on this by using larger A/l at high L, to pro- 
duce more loops there (using z/l automatically 
renormalizes their contributions). We work in 
D — 3(4:) and use an {uj{Xf_i — I)} ensemble of up 
to 600000(250000) loops on lattices oiN^NyN^Nt 
with N^, = 12, 16, 24 and 32, p.b.c. See Fig. |l|. 




Figure 1. RW vl (§) vs L in = 3 (diamonds) 
and 4 (crosses). Straight lines represent L^^/^. 

Special configurations For a test we use: 
a) Homogeneous field In continuum [0: 



£(6+,6_) = 
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(with subtraction at T = 0). With 

b^ = 2Trk^/N^, k^, : integer, (5) 
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, l<n^<N^ (6) 



with = (all other links 1) realizes a con- 
stant, Abelian magnetic field b ~ b^ — bx + 
by, b- = 0. Alternatively, with iV^ = A^^ = A^, 
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(all other links 1) realizes an SU(2) configuration 
of constant topological charge density q = 
(here ky^t = k; b± = b = bx = h). We calculate 
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From (|^,^) this should be 

F^^ir.{L,b) ^u/sinhu, u = bL/2D, (9) 
Ftop.(i,fo) = uVsinh^u, u^bL/AD (10) 

for the first (||) and second (Q) case, respectively. 
For small fields the points scale well with the field 
and lattice size and reproduce the continuum re- 
sults (^ or ([l^). For large fields lattice artifacts 
are seen on small lattices. See Fig. ||. 
h) Magnetic field in half space Consider a mag- 
netic field (^ in _D = 3 pointing into z-direction 
in half space I < x < = ^. Then: 

C{Li,L2,b,m,x) ^ -e^'^VLF{L,b,x), 



1 

F(L,5,x) = — Vtr W Ui 



(11) 
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The lattice is 64 x 32^, m > 0, uf'ix) are centered 
at X. On physical grounds we expect for L |^ 



/(0,oo,^,e) cxe-^^(^)« for ^>0 
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Here a;o = + 0.5 (^ — 0) is the separation be- 
tween regions. We have systematic errors due to 
< Ai,2 < c» and discretization artifacts. The 
slope M appears reasonably stable against these 
uncertainties, however. See Fig. ||. We also ob- 
serve scaling (for small b and m) and smooth de- 
pendence on m. We obtain M(0) — 1.7 ± 0.2, in 
apparent disagreement with n. See B. 




Figure 2. F{L, b) vs x = bL in D ^ 4, for U(l) 
(|) (6 in units of Air/N, upper plot) and SU(2) 
case (6 in units of 2tt/N, lower plot). The 
lines represent the continuum results (||), ([lO|). 



Discussion The above results indicate that 
using statistical ensembles of RW loops allows one 
to calculate the matter free energy in gauge field 
backgrounds where large loops may be relevant. 
This method could be applied to non-trivial con- 
figurations such as instantons, sphalerons, bounce 
solutions. It also could be used to calculate the 
quenched free energy in a quantum vacuum [ pO[ . 
Non-zero spin can also be handled by adding a 
spin term in the worldline Lagrangian containing 
the outer field Q (for fermions it is advantageous 
to use second order formalism - suggestion by W. 
Wetzel). See also Of course one can use di- 
rect numerical methods for determinants. The 
above method may, however, allow more insight 
into the relevant phenomena and scales involved. 
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